Evolutionary analysis reveals the role of a non-catalytic domain of peptidyl arginine deiminase 2 in transcriptional regulation.

Peptidyl arginine deiminases (PADIs) catalyze protein citrullination, a post-translational conversion of arginine to citrulline. The most widely expressed member of this family, PADI2, regulates cellular processes that impact several diseases. We hypothesized that we could gain new insights into PADI2 function through a systematic evolutionary and structural analysis. Here, we identify 20 positively selected PADI2 residues, 16 of which are structurally exposed and maintain PADI2 interactions with cognate proteins. Many of these selected residues reside in non-catalytic regions of PADI2. We validate the importance of a prominent loop in the middle domain that encompasses PADI2 L162, a residue under positive selection. This site is essential for interaction with the transcription elongation factor (P-TEFb) and mediates the active transcription of the oncogenes c-MYC, and CCNB1, as well as impacting cellular proliferation. These insights could be key to understanding and addressing the role of the PADI2 c-MYC axis in cancer progression.

Shared Proteins and Pathways of Cardiovascular and Cognitive Diseases: Relation to Vascular Cognitive Impairment.

One of the primary goals of systems medicine is the detection of putative proteins and pathways involved in disease progression and pathological phenotypes. Vascular cognitive impairment (VCI) is a heterogeneous condition manifesting as cognitive impairment resulting from vascular factors. The precise mechanisms underlying this relationship remain unclear, which poses challenges for experimental research. Here, we applied computational approaches like systems biology to unveil and select relevant proteins and pathways related to VCI by studying the crosstalk between cardiovascular and cognitive diseases. In addition, we specifically included signals related to oxidative stress, a common etiologic factor tightly linked to aging, a major determinant of VCI. Our results show that pathways associated with oxidative stress are quite relevant, as most of the prioritized vascular cognitive genes and proteins were enriched in these pathways. Our analysis provided a short list of proteins that could be contributing to VCI: DOLK, TSC1, ATP1A1, MAPK14, YWHAZ, CREB3, HSPB1, PRDX6, and LMNA. Moreover, our experimental results suggest a high implication of glycative stress, generating oxidative processes and post-translational protein modifications through advanced glycation end-products (AGEs). We propose that these products interact with their specific receptors (RAGE) and Notch signaling to contribute to the etiology of VCI.

A quantitative systems pharmacology model for certolizumab pegol treatment in moderate-to-severe psoriasis.

Background: Psoriasis is a chronic immune-mediated inflammatory systemic disease with skin manifestations characterized by erythematous, scaly, itchy and/or painful plaques resulting from hyperproliferation of keratinocytes. Certolizumab pegol [CZP], a PEGylated antigen binding fragment of a humanized monoclonal antibody against TNF-alpha, is approved for the treatment of moderate-to-severe plaque psoriasis. Patients with psoriasis present clinical and molecular variability, affecting response to treatment. Herein, we utilized an in silico approach to model the effects of CZP in a virtual population (vPop) with moderate-to-severe psoriasis. Our proof-of-concept study aims to assess the performance of our model in generating a vPop and defining CZP response variability based on patient profiles. Methods: We built a quantitative systems pharmacology (QSP) model of a clinical trial-like vPop with moderate-to-severe psoriasis treated with two dosing schemes of CZP (200 mg and 400 mg, both every two weeks for 16 weeks, starting with a loading dose of CZP 400 mg at weeks 0, 2, and 4). We applied different modelling approaches: (i) an algorithm to generate vPop according to reference population values and comorbidity frequencies in real-world populations; (ii) physiologically based pharmacokinetic (PBPK) models of CZP dosing schemes in each virtual patient; and (iii) systems biology-based models of the mechanism of action (MoA) of the drug. Results: The combination of our different modelling approaches yielded a vPop distribution and a PBPK model that aligned with existing literature. Our systems biology and QSP models reproduced known biological and clinical activity, presenting outcomes correlating with clinical efficacy measures. We identified distinct clusters of virtual patients based on their psoriasis-related protein predicted activity when treated with CZP, which could help unravel differences in drug efficacy in diverse subpopulations. Moreover, our models revealed clusters of MoA solutions irrespective of the dosing regimen employed. Conclusion: Our study provided patient specific QSP models that reproduced clinical and molecular efficacy features, supporting the use of computational methods as modelling strategy to explore drug response variability. This might shed light on the differences in drug efficacy in diverse subpopulations, especially useful in complex diseases such as psoriasis, through the generation of mechanistically based hypotheses.

Schlafen 12 restricts HIV-1 latency reversal by a codon-usage dependent post-transcriptional block in CD4+ T cells.

Latency is a major barrier towards virus elimination in HIV-1-infected individuals. Yet, the mechanisms that contribute to the maintenance of HIV-1 latency are incompletely understood. Here we describe the Schlafen 12 protein (SLFN12) as an HIV-1 restriction factor that establishes a post-transcriptional block in HIV-1-infected cells and thereby inhibits HIV-1 replication and virus reactivation from latently infected cells. The inhibitory activity is dependent on the HIV-1 codon usage and on the SLFN12 RNase active sites. Within HIV-1-infected individuals, SLFN12 expression in PBMCs correlated with HIV-1 plasma viral loads and proviral loads suggesting a link with the general activation of the immune system. Using an RNA FISH-Flow HIV-1 reactivation assay, we demonstrate that SLFN12 expression is enriched in infected cells positive for HIV-1 transcripts but negative for HIV-1 proteins. Thus, codon-usage dependent translation inhibition of HIV-1 proteins participates in HIV-1 latency and can restrict the amount of virus release after latency reversal.

The P4-ATPase Drs2 interacts with and stabilizes the multisubunit tethering complex TRAPPIII in yeast.

Multisubunit Tethering Complexes (MTCs) are a set of conserved protein complexes that tether vesicles at the acceptor membrane. Interactions with other components of the trafficking machinery regulate MTCs through mechanisms that are partially understood. Here, we systematically investigate the interactome that regulates MTCs. We report that P4-ATPases, a family of lipid flippases, interact with MTCs that participate in the anterograde and retrograde transport at the Golgi, such as TRAPPIII. We use the P4-ATPase Drs2 as a paradigm to investigate the mechanism and biological relevance of this interplay during transport of Atg9 vesicles. Binding of Trs85, the sole-specific subunit of TRAPPIII, to the N-terminal tail of Drs2 stabilizes TRAPPIII on membranes loaded with Atg9 and is required for Atg9 delivery during selective autophagy, a role that is independent of P4-ATPase canonical functions. This mechanism requires a conserved I(S/R)TTK motif that also mediates the interaction of the P4-ATPases Dnf1 and Dnf2 with MTCs, suggesting a broader role of P4-ATPases in MTC regulation.

A Genome-Wide Functional Screen Identifies Enhancer and Protective Genes for Amyloid Beta-Peptide Toxicity.

Alzheimer's disease (AD) is known to be caused by amyloid ß-peptide (Aß) misfolded into ß-sheets, but this knowledge has not yet led to treatments to prevent AD. To identify novel molecular players in Aß toxicity, we carried out a genome-wide screen in Saccharomyces cerevisiae, using a library of 5154 gene knock-out strains expressing Aß1-42. We identified 81 mammalian orthologue genes that enhance Aß1-42 toxicity, while 157 were protective. Next, we performed interactome and text-mining studies to increase the number of genes and to identify the main cellular functions affected by Aß oligomers (oAß). We found that the most affected cellular functions were calcium regulation, protein translation and mitochondrial activity. We focused on SURF4, a protein that regulates the store-operated calcium channel (SOCE). An in vitro analysis using human neuroblastoma cells showed that SURF4 silencing induced higher intracellular calcium levels, while its overexpression decreased calcium entry. Furthermore, SURF4 silencing produced a significant reduction in cell death when cells were challenged with oAß1-42, whereas SURF4 overexpression induced Aß1-42 cytotoxicity. In summary, we identified new enhancer and protective activities for Aß toxicity and showed that SURF4 contributes to oAß1-42 neurotoxicity by decreasing SOCE activity.

Amyloid ß-Peptide Causes the Permanent Activation of CaMKIIα through Its Oxidation.

Alzheimer's disease (AD) is characterised by the presence of extracellular amyloid plaques in the brain. They are composed of aggregated amyloid beta-peptide (Aß) misfolded into beta-sheets which are the cause of the AD memory impairment and dementia. Memory depends on the hippocampal formation and maintenance of synapses by long-term potentiation (LTP), whose main steps are the activation of NMDA receptors, the phosphorylation of CaMKIIα and the nuclear translocation of the transcription factor CREB. It is known that Aß oligomers (oAß) induce synaptic loss and impair the formation of new synapses. Here, we have studied the effects of oAß on CaMKIIα. We found that oAß produce reactive oxygen species (ROS), that induce CaMKIIα oxidation in human neuroblastoma cells as we assayed by western blot and immunofluorescence. Moreover, this oxidized isoform is significantly present in brain samples from AD patients. We found that the oxidized CaMKIIα is active independently of the binding to calcium/calmodulin, and that CaMKIIα phosphorylation is mutually exclusive with CaMKIIα oxidation as revealed by immunoprecipitation and western blot. An in silico modelling of the enzyme was also performed to demonstrate that oxidation induces an activated state of CaMKIIα. In brains from AD transgenic models of mice and in primary cultures of murine hippocampal neurons, we demonstrated that the oxidation of CaMKIIα induces the phosphorylation of CREB and its translocation to the nucleus to promote the transcription of ARC and BDNF. Our data suggests that CaMKIIα oxidation would be a pro-survival mechanism that is triggered when a noxious stimulus challenges neurons as do oAß.

Differential regulation of insulin signalling by monomeric and oligomeric amyloid beta-peptide.

Alzheimer's disease and Type 2 diabetes are pathological processes associated to ageing. Moreover, there are evidences supporting a mechanistic link between Alzheimer's disease and insulin resistance (one of the first hallmarks of Type 2 diabetes). Regarding Alzheimer's disease, amyloid ß-peptide aggregation into ß-sheets is the main hallmark of Alzheimer's disease. At monomeric state, amyloid ß-peptide is not toxic but its function in brain, if any, is unknown. Here we show, by in silico study, that monomeric amyloid ß-peptide 1-40 shares the tertiary structure with insulin and is thereby able to bind and activate insulin receptor. We validated this prediction experimentally by treating human neuroblastoma cells with increasing concentrations of monomeric amyloid ß-peptide 1-40. Our results confirm that monomeric amyloid ß-peptide 1-40 activates insulin receptor autophosphorylation, triggering downstream enzyme phosphorylations and the glucose Transporter 4 translocation to the membrane. On the other hand, neuronal insulin resistance is known to be associated to Alzheimer's disease since early stages. We thus modelled the docking of oligomeric amyloid ß-peptide 1-40 to insulin receptor. We found that oligomeric amyloid ß-peptide 1-40 blocks insulin receptor, impairing its activation. It was confirmed in vitro by observing the lack of insulin receptor autophosphorylation, and also the impairment of insulin-induced intracellular enzyme activations and the glucose Transporter 4 translocation to the membrane. By biological system analysis, we have carried out a mathematical model recapitulating the process that turns amyloid ß-peptide binding to insulin receptor from the physiological to the pathophysiological regime. Our results suggest that monomeric amyloid ß-peptide 1-40 contributes to mimic insulin effects in the brain, which could be good when neurons have an extra requirement of energy beside the well-known protective effects on insulin intracellular signalling, while its accumulation and subsequent oligomerization blocks the insulin receptor producing insulin resistance and compromising neuronal metabolism and protective pathways.

Melanoma RBPome identification reveals PDIA6 as an unconventional RNA-binding protein involved in metastasis.

RNA-binding proteins (RBPs) have been relatively overlooked in cancer research despite their contribution to virtually every cancer hallmark. Here, we use RNA interactome capture (RIC) to characterize the melanoma RBPome and uncover novel RBPs involved in melanoma progression. Comparison of RIC profiles of a non-tumoral versus a metastatic cell line revealed prevalent changes in RNA-binding capacities that were not associated with changes in RBP levels. Extensive functional validation of a selected group of 24 RBPs using five different in vitro assays unveiled unanticipated roles of RBPs in melanoma malignancy. As proof-of-principle we focused on PDIA6, an ER-lumen chaperone that displayed a novel RNA-binding activity. We show that PDIA6 is involved in metastatic progression, map its RNA-binding domain, and find that RNA binding is required for PDIA6 tumorigenic properties. These results exemplify how RIC technologies can be harnessed to uncover novel vulnerabilities of cancer cells.

Characterization of Novel Pathogenic Variants Leading to Caspase-8 Cleavage-Resistant RIPK1-Induced Autoinflammatory Syndrome.

Pathogenic RIPK1 variants have been described as the cause of two different inborn errors of immunity. Biallelic loss-of-function variants cause the recessively inherited RIPK1 deficiency, while monoallelic variants impairing the caspase-8-mediated RIPK1 cleavage provoke a novel autoinflammatory disease (AID) called cleavage-resistant RIPK1-induced autoinflammatory (CRIA) syndrome. The aim of this study was to characterize the pathogenicity of two novel RIPK1 variants located at the cleavage site of caspase-8 detected in patients with dominantly-inherited, early-onset undefined AID. RIPK1 genotyping was performed by Sanger and next-generation sequencing. Clinical and analytical data were collected from medical charts, and in silico and in vitro assays were performed to evaluate the functional consequences. Genetic analyses identified two novel heterozygous RIPK1 variants at the caspase-8 cleavage site (p.Leu321Arg and p.Asp324Gly), which displayed a perfect intrafamilial phenotype-genotype segregation following a dominant inheritance pattern. Structural analyses suggested that these variants disrupt the normal RIPK1 structure, probably making it less accessible to and/or less cleavable by caspase-8. In vitro experiments confirmed that the p.Leu321Arg and p.Asp324Gly RIPK1 variants were resistant to caspase-8-mediated cleavage and induced a constitutive activation of necroptotic pathway in a similar manner that previously characterized RIPK1 variants causing CRIA syndrome. All these results strongly supported the pathogenicity of the two novel RIPK1 variants and the diagnosis of CRIA syndrome in all enrolled patients. Moreover, the evidences here collected expand the phenotypic and genetic diversity of this recently described AID, and provide interesting data about effectiveness of treatments that may benefit future patients.

First Description of Late-Onset Autoinflammatory Disease Due to Somatic NLRC4 Mosaicism.

OBJECTIVE: Autoinflammatory diseases are inherited disorders of innate immunity that usually start during childhood. However, several recent reports have described an increasing number of patients with autoinflammatory disease starting in adulthood. This study was undertaken to investigate the underlying cause of a case of late-onset uncharacterized autoinflammatory disease. METHODS: Genetics studies were performed using Sanger sequencing and next-generation sequencing (NGS) methods. In silico, in vitro, and ex vivo analyses were performed to determine the functional consequences of the detected variant. RESULTS: We studied a 57-year-old woman who at the age of 47 years began to have recurrent episodes of fever, myalgias, arthralgias, diffuse abdominal pain, diarrhea, adenopathies, and systemic inflammation, which were relatively well controlled with anti-interleukin-1 (anti-IL-1) drugs. NGS analyses did not detect germline variants in any of the known autoinflammatory disease-associated genes, but they identified the p.Ser171Phe NLRC4 variant in unfractionated blood, with an allele fraction (2-4%) compatible with gene mosaicism. Structural modeling analyses suggested that this missense variant might favor the open, active conformation of the NLRC4 protein, and in vitro and ex vivo analyses confirmed its propensity to oligomerize and activate the NLRC4 inflammasome, with subsequent overproduction of IL-18. CONCLUSION: Our findings indicate that the postzygotic p.Ser171Phe NLRC4 variant is a plausible cause of the disease in the enrolled patient. Functional and structural studies clearly support, for the first time, its gain-of-function behavior, consistent with previously reported NLRC4 pathogenic variants. These novel findings should be considered in the diagnostic evaluation of patients with adult-onset uncharacterized autoinflammatory disease.

Find and cut-and-transfer (FiCAT) mammalian genome engineering.

While multiple technologies for small allele genome editing exist, robust technologies for targeted integration of large DNA fragments in mammalian genomes are still missing. Here we develop a gene delivery tool (FiCAT) combining the precision of a CRISPR-Cas9 (find module), and the payload transfer efficiency of an engineered piggyBac transposase (cut-and-transfer module). FiCAT combines the functionality of Cas9 DNA scanning and targeting DNA, with piggyBac donor DNA processing and transfer capacity. PiggyBac functional domains are engineered providing increased on-target integration while reducing off-target events. We demonstrate efficient delivery and programmable insertion of small and large payloads in cellulo (human (Hek293T, K-562) and mouse (C2C12)) and in vivo in mouse liver. Finally, we evolve more efficient versions of FiCAT by generating a targeted diversity of 394,000 variants and undergoing 4 rounds of evolution. In this work, we develop a precise and efficient targeted insertion of multi kilobase DNA fragments in mammalian genomes.

Methods to Develop an in silico Clinical Trial: Computational Head-to-Head Comparison of Lisdexamfetamine and Methylphenidate.

Regulatory agencies encourage computer modeling and simulation to reduce the time and cost of clinical trials. Although still not classified in formal guidelines, system biology-based models represent a powerful tool for generating hypotheses with great molecular detail. Herein, we have applied a mechanistic head-to-head in silico clinical trial (ISCT) between two treatments for attention-deficit/hyperactivity disorder, to wit lisdexamfetamine (LDX) and methylphenidate (MPH). The ISCT was generated through three phases comprising (i) the molecular characterization of drugs and pathologies, (ii) the generation of adult and children virtual populations (vPOPs) totaling 2,600 individuals and the creation of physiologically based pharmacokinetic (PBPK) and quantitative systems pharmacology (QSP) models, and (iii) data analysis with artificial intelligence methods. The characteristics of our vPOPs were in close agreement with real reference populations extracted from clinical trials, as did our PBPK models with in vivo parameters. The mechanisms of action of LDX and MPH were obtained from QSP models combining PBPK modeling of dosing schemes and systems biology-based modeling technology, i.e., therapeutic performance mapping system. The step-by-step process described here to undertake a head-to-head ISCT would allow obtaining mechanistic conclusions that could be extrapolated or used for predictions to a certain extent at the clinical level. Altogether, these computational techniques are proven an excellent tool for hypothesis-generation and would help reach a personalized medicine.

Sensing low intracellular potassium by NLRP3 results in a stable open structure that promotes inflammasome activation.

The NLRP3 inflammasome is activated by a wide range of stimuli and drives diverse inflammatory diseases. The decrease of intracellular K+ concentration is a minimal upstream signal to most of the NLRP3 activation models. Here, we found that cellular K+ efflux induces a stable structural change in the inactive NLRP3, promoting an open conformation as a step preceding activation. This conformational change is facilitated by the specific NLRP3 FISNA domain and a unique flexible linker sequence between the PYD and FISNA domains. This linker also facilitates the ensemble of NLRP3PYD into a seed structure for ASC oligomerization. The introduction of the NLRP3 PYD-linker-FISNA sequence into NLRP6 resulted in a chimeric receptor able to be activated by K+ efflux­specific NLRP3 activators and promoted an in vivo inflammatory response to uric acid crystals. Our results establish that the amino-terminal sequence between PYD and NACHT domain of NLRP3 is key for inflammasome activation.

Comparative Analysis of Mammal Genomes Unveils Key Genomic Variability for Human Life Span.

The enormous mammal's lifespan variation is the result of each species' adaptations to their own biological trade-offs and ecological conditions. Comparative genomics have demonstrated that genomic factors underlying both, species lifespans and longevity of individuals, are in part shared across the tree of life. Here, we compared protein-coding regions across the mammalian phylogeny to detect individual amino acid (AA) changes shared by the most long-lived mammals and genes whose rates of protein evolution correlate with longevity. We discovered a total of 2,737 AA in 2,004 genes that distinguish long- and short-lived mammals, significantly more than expected by chance (P = 0.003). These genes belong to pathways involved in regulating lifespan, such as inflammatory response and hemostasis. Among them, a total 1,157 AA showed a significant association with maximum lifespan in a phylogenetic test. Interestingly, most of the detected AA positions do not vary in extant human populations (81.2%) or have allele frequencies below 1% (99.78%). Consequently, almost none of these putatively important variants could have been detected by genome-wide association studies. Additionally, we identified four more genes whose rate of protein evolution correlated with longevity in mammals. Crucially, SNPs located in the detected genes explain a larger fraction of human lifespan heritability than expected, successfully demonstrating for the first time that comparative genomics can be used to enhance interpretation of human genome-wide association studies. Finally, we show that the human longevity-associated proteins are significantly more stable than the orthologous proteins from short-lived mammals, strongly suggesting that general protein stability is linked to increased lifespan.

Head to head evaluation of second generation ALK inhibitors brigatinib and alectinib as first-line treatment for ALK+ NSCLC using an in silico systems biology-based approach.

Around 3-7% of patients with non-small cell lung cancer (NSCLC), which represent 85% of diagnosed lung cancers, have a rearrangement in the ALK gene that produces an abnormal activity of the ALK protein cell signaling pathway. The developed ALK tyrosine kinase inhibitors (TKIs), such as crizotinib, ceritinib, alectinib, brigatinib and lorlatinb present good performance treating ALK+ NSCLC, although all patients invariably develop resistance due to ALK secondary mutations or bypass mechanisms. In the present study, we compare the potential differences between brigatinib and alectinib's mechanisms of action as first-line treatment for ALK+ NSCLC in a systems biology-based in silico setting. Therapeutic performance mapping system (TPMS) technology was used to characterize the mechanisms of action of brigatinib and alectinib and the impact of potential resistances and drug interferences with concomitant treatments. The analyses indicate that brigatinib and alectinib affect cell growth, apoptosis and immune evasion through ALK inhibition. However, brigatinib seems to achieve a more diverse downstream effect due to a broader cancer-related kinase target spectrum. Brigatinib also shows a robust effect over invasiveness and central nervous system metastasis-related mechanisms, whereas alectinib seems to have a greater impact on the immune evasion mechanism. Based on this in silico head to head study, we conclude that brigatinib shows a predicted efficacy similar to alectinib and could be a good candidate in a first-line setting against ALK+ NSCLC. Future investigation involving clinical studies will be needed to confirm these findings. These in silico systems biology-based models could be applied for exploring other unanswered questions.

Human Albumin Impairs Amyloid ß-peptide Fibrillation Through its C-terminus: From docking Modeling to Protection Against Neurotoxicity in Alzheimer's disease.

Alzheimer's disease (AD) is a neurodegenerative process characterized by the accumulation of extracellular deposits of amyloid ß-peptide (Aß), which induces neuronal death. Monomeric Aß is not toxic but tends to aggregate into ß-sheets that are neurotoxic. Therefore to prevent or delay AD onset and progression one of the main therapeutic approaches would be to impair Aß assembly into oligomers and fibrils and to promote disaggregation of the preformed aggregate. Albumin is the most abundant protein in the cerebrospinal fluid and it was reported to bind Aß impeding its aggregation. In a previous work we identified a 35-residue sequence of clusterin, a well-known protein that binds Aß, that is highly similar to the C-terminus (CTerm) of albumin. In this work, the docking experiments show that the average binding free energy of the CTerm-Aß1-42 simulations was significantly lower than that of the clusterin-Aß1-42 binding, highlighting the possibility that the CTerm retains albumin's binding properties. To validate this observation, we performed in vitro structural analysis of soluble and aggregated 1 µM Aß1-42 incubated with 5 µM CTerm, equimolar to the albumin concentration in the CSF. Reversed-phase chromatography and electron microscopy analysis demonstrated a reduction of Aß1-42 aggregates when the CTerm was present. Furthermore, we treated a human neuroblastoma cell line with soluble and aggregated Aß1-42 incubated with CTerm obtaining a significant protection against Aß-induced neurotoxicity. These in silico and in vitro data suggest that the albumin CTerm is able to impair Aß aggregation and to promote disassemble of Aß aggregates protecting neurons.

The exonuclease Xrn1 activates transcription and translation of mRNAs encoding membrane proteins.

The highly conserved 5'-3' exonuclease Xrn1 regulates gene expression in eukaryotes by coupling nuclear DNA transcription to cytosolic mRNA decay. By integrating transcriptome-wide analyses of translation with biochemical and functional studies, we demonstrate an unanticipated regulatory role of Xrn1 in protein synthesis. Xrn1 promotes translation of a specific group of transcripts encoding membrane proteins. Xrn1-dependence for translation is linked to poor structural RNA contexts for translation initiation, is mediated by interactions with components of the translation initiation machinery and correlates with an Xrn1-dependence for mRNA localization at the endoplasmic reticulum, the translation compartment of membrane proteins. Importantly, for this group of mRNAs, Xrn1 stimulates transcription, mRNA translation and decay. Our results uncover a crosstalk between the three major stages of gene expression coordinated by Xrn1 to maintain appropriate levels of membrane proteins.

Arginine Citrullination at the C-Terminal Domain Controls RNA Polymerase II Transcription.

The post-translational modification of key residues at the C-terminal domain of RNA polymerase II (RNAP2-CTD) coordinates transcription, splicing, and RNA processing by modulating its capacity to act as a landing platform for a variety of protein complexes. Here, we identify a new modification at the CTD, the deimination of arginine and its conversion to citrulline by peptidyl arginine deiminase 2 (PADI2), an enzyme that has been associated with several diseases, including cancer. We show that, among PADI family members, only PADI2 citrullinates R1810 (Cit1810) at repeat 31 of the CTD. Depletion of PADI2 or loss of R1810 results in accumulation of RNAP2 at transcription start sites, reduced gene expression, and inhibition of cell proliferation. Cit1810 is needed for interaction with the P-TEFb (positive transcription elongation factor b) kinase complex and for its recruitment to chromatin. In this way, CTD-Cit1810 favors RNAP2 pause release and efficient transcription in breast cancer cells.